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Nonequilibrium dynamics of the ±J Ising, the XY, and the Heisenberg spin-glass models are in- 
vestigated in three dimensions. A nonequilibrium dynamic exponent is calculated from the dynamic 
correlation length. The spin-glass dynamic exponent continuously depends on the temperature. 
There is no anomaly at the critical temperature as is recently reported by Katzgraber and Campbell. 
On the other hand, the chiral-glass dynamic exponent takes a constant value above the spin-glass 
transition temperature (Tsg), and becomes temperature-dependent below Tsg. The finite-time scal- 
ing analyses on the spin- and the chiral-glass susceptibility are performed using the temperature 
dependence of the dynamic exponent. The critical temperature and the critical exponents are es- 
timated with high precisions. In the Ising model, Tsg = 1.18(2), Usg — 1.75(10), 7sg = 3.4(2), and 
z,g = 5.1(1). In the model, T^g = 0.435(15), u^g = 1.25(15), 7sg = 2.30(25), and = 5.0(1) 
for the spin-glass transition, and Teg — 0.445(15), Vcg = 1.05(10), 7cg = 1.65(20), and Zcg = 5.1(2) 
for the chiral-glass transition. In the Heisenberg model, Tsg = 0.18(1), Vsg = 1.55(15), 7sg — 2.9(2), 
and Zsg = 5.1(1) for the spin-glass transition, and Teg — 0.23(1), Ucg — 0.70(15), 7cg — 0.65(15), and 
Zcg = 5.0(2) for the chiral-glass transition. A difference of the spin- and the chiral-glass transition 
temperatures is resolved in the Heisenberg model. The dynamic critical exponent takes almost the 
same value for all transitions. It suggests that the spin-glass and the chiral-glass transitions in three 
dimensions are dynamically universal. 



I. INTRODUCTION 

Randomness and frustration produce nontrivial phe- 
nomena in spin glassesiiSi^. There is no spatial or- 
der in the spin-glass phase. The spins are ordered 
only in the time direction: spins randomly freeze be- 
low the critical temperature. Growth of the spin- 
glass correlations is very slow. Most of interesting 
observations in experiments and in computer simula- 
tions are essentially nonequilibrium phenomena. Ag- 
ing, memory and rejuvenation effects are the typical 
examples 4i^*S4i It has become clear that a domain size of 

frozen spins plays an important role in explaining these 
phenomenai«i^i"i"ii^ii3i^4ii5ii'^-i^iiS-i9,2Qi2i,22,23 q^^^^^^ 

of time evolution of the domain size may determine the 
dynamics of spin glasses. 

In the second-order phase transitions, a correlation 
length ^ and a correlation time t diverge at the criti- 
cal temperature. The dynamic scaling relation holds as 
T ~ .^^ with the dynamic critical exponent z. If we ob- 
serve the correlation length as a function of time, it also 
diverges algebraically with the same exponent. It is writ- 
ten as S,{t) ^ t^^^. Consequently, the magnetic suscepti- 
bility xit) diverges with an exponent 'y / zv. 

The dynamic scaling relation also holds in spin glasses 
at the critical temperature. Huse^i remarked that a time 
evolution of the spin-glass susceptibility Xsg(i) diverges 
algebraically with an exponent j/zv in the Ising spin- 
glass model. Blundell et alm^ found that a relaxation 
function of the Binder parameter gsg also diverges with 
an exponent 1/z. The spin-glass correlation length is also 
considered to diverge algebraically at the critical temper- 
ature. 

Below the spin-glass transition temperature, Tgg, the 



algebraic divergence of the spin-glass correlation length 

is considered to continue in the nonequilibrium process. 

An exponent of the divergence linearly depends on the 
temperaturei^^ii5,^6,i7,i8,i9 
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Here, we call Zs{T) a spin-glass nonequilibrium dynamic 
exponent. It is defined only in the nonequilibrium relax- 
ation process. It equals to the dynamic critical exponent 
of the spin glass transition, Zgg, at the critical tempera- 
ture. 

Recently, Katzgraber and CampbellS^ remarked that 
the relation continues to hold above Tsg, and that 
there is no anomaly at the critical temperature. The 
model they investigated is the Gaussian-bond Ising spin- 
glass model and the gauge-glass model in three and four 
dimensions. It is very important whether or not this 
finding is common to various spin-glass models. In regu- 
lar systems, the dynamic critical exponent z governs the 
algebraic growth of the correlation length in the nonequi- 
librium relaxation process above the critical temperature. 
The temperature dependence of z above Tsg is out of our 
assumption regarding the dynamics of spin glasses. 

The dynamic scaling relation in the nonequilibrium re- 
laxation process of the spin-glass models suggests that 
we can study the spin-glass phase transition using the 
nonequilibrium relaxation mcthodi2SiSL22iSS It is a pow- 
erful tool to estimate the critical temperature and the 
critical exponent with high precisions. It has been ap- 
plied to slow dynamic systems with frustration and/or 
randomness i^^i^^i^°i^^i'^^i^^i^"^i^^ A key concept of the 
method is an exchange of time and size by the dynamic 
scaling relation, t ^ L^. Here, a single exponent z is 
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considered to satisfy the relation above the critical tem- 
perature. This assumption may be broken by the relation 
Therefore, we must reconsider the nonequilibrium 
relaxation analysis in the spin-glass problems taking into 
account the relation Then, the finite-time-scaling 

results may be influenced. 

In this paper we carry out extensive Monte Carlo sim- 
ulations. It is made clear that the relation ^ is satisfied 
for the spin-glass transitions in all models. It contin- 
ues regardless the critical temperature without anomaly. 
On the other hand, the temperature dependence of the 
nonequilibrium dynamic exponent of the chiral-glass cor- 
relation length only appears below Tsg- It takes a con- 
stant value above Tgg as in the regular systems. The 
finite-time scaling analyses are performed using the tem- 
perature dependence of the dynamic exponent. The tran- 
sition temperature and the critical exponents are esti- 
mated with higher precisions compared to our previous 
results. We give a conclusion to a problem of whether or 
not the spin- and the chiral-glass transitions occur simul- 
taneously in the XY and the Heisenberg models. 

The present paper is organized as follows. Section UTI 
describes a model and physical quantities we observe in 
simulations. Section IIIII explains our numerical method 
and some technical details. Results of the Ising, the XY, 
and the Heisenberg models are presented in Sec. IIVI Ivl 
and I VII respectively. Section Ivlll is devoted to summary 
and discussions. 



II. MODEL AND OBSERVABLES 
A. Model 

We study the ±J spin-glass model in a simple cubic 
lattice. 
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The sum runs over all the nearest-neighbor spin pairs 
(ij). The interactions, J^, take the two values +J and 
—J with the same probability. The temperature, T, is 
scaled by J. Linear lattice size is denoted by L. A total 
number of spins is N = L x L X {L + 1), and skewed pe- 
riodic boimdary conditions are imposed. The spins are 
updated by the single-spin-flip algorithm using the con- 
ventional Metropolis probability. We treat three types of 
spin dimensions: the Ising model with Si = Siz, the XY 
model with Si = iSix,Siy), and the Heisenberg model 
with Si ~ {Six, Siy, Siz). In the model, each spin is 
defined by an angle in the Xy plane: Si ~ {cos 9i, sin 6';). 
The angle 6i takes continuous values in [0, 27r); however, 
we divide the interval into 1024 discrete states in this 
study. The effect of the discretization is checked by com- 
paring data of 1024-state simulations with those of 2048- 
state simulations. It is found to be negligible. 

The physical quantities that are calculated in our sim- 
ulations arc the spin-glass susceptibility, ^sg, the chiral- 



glass susceptibility, Xcg, and the spin- and the chiral-glass 
correlation functions, fsg and /eg, from which we estimate 
the correlation lengths, ^sg and ^cg- 



B. Spin-glass susceptibility 

The spin-glass susceptibility is defined by the following 
expression. 



1 

Xsg = ^ 
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The thermal average is denoted by {■■■), and the 
random-bond configurational average is denoted by 
[•••]c. The thermal average is replaced by an average 
over independent real replicas that consist of different 
thermal ensembles: 
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A replica number is denoted by m. The superscripts A 
is the replica index. We prepare m real replicas for each 
random bond with different initial spin configurations. 
Each replica is updated using a different random num- 
ber sequence. The thermal average is taken only by this 
replica average in our numerical scheme. A replica num- 
ber controls an accuracy of the thermal average. 

An overlap between two replicas, A and B, is defined 

by 



1 V 9(^4) q(-S) 
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Here, subscripts and v represent three spin directions 
and z. The spin-glass susceptibility is estimated 
using the overlap function: 
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Similarly, the chiral-glass susceptibility in the XY 
mode l?^i?T-?? is defined by 
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mym — 1) 
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where 
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A local chirality variable, the z component of the 

vector chirality. It is defined at each square plaquctte. 
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a, that consists of the nearest-neighbor bonds, (jk). In 
the Heisenberg model, we use the scalar chiralitji22il2i^, 
which is defined by 



(A) 



(10) 



where denotes a unit lattice vector along the /i axis. 
The overlap function and the chiral-glass susceptibility 
are defined in the same manner as in the XY model. 



C. Dynamic correlation length 

A spin-glass correlation function is defined by the fol- 
lowing expression 1)11(1 . 
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Here, q^^ii) - ^.^^ 



5*1^^ 5-^^ is the spin overlap on the ith 
site. When a replica number is two, it is equivalent to the 
four-point correlation function as is shown in Eq. (|12|1 . 
Since m is large in this study, it is very time-consuming 
to take an average over m(m — l)/2 different overlap 
functions. Therefore, we use another expression H13|) for 
/sg(r). For a given distance r and a site i, we calculate 
a spin correlation function for replica A, and store it in 
a memory array. Then, the replica average is taken and 
the obtained value is squared. Taking a summation of 
the squared value over lattice sites i we obtain the corre- 
lation function fsgir) for a given distance r. Changing a 
value of r we do the same procedure to evaluate /sg (r) . 
Total amount of calculation is much smaller in this pro- 
cedure because the maximum value of r = _L/2 is much 
smaller than m{m — l)/2. 

We obtain a chiral-glass correlation function in the 
same manner replacing the local spin variables with the 
local chirality variables. Here, we only consider correla- 
tions between two parallel square plaquettes in the XY 
model. A pair of three neighboring spins in the same 
direction is considered in the Heisenberg model. We es- 
timate the spin- and the chiral-glass correlation lengths, 
^sg and ^cgj from the correlation functions using the fol- 
lowing expression. 



fsg/cg{r) = /oexp ^- 



Csg/c 



(14) 



A problem in this expression is a prefactor /q. We 
may suppose several forms for /g. For example, /o ~ 
^2-d-T]_ However, the obtained value may depend on a 



value of r], which is not known yet. Therefore, we take 
a rather brute-force approach as explained below. Then, 
the relaxation functions of the correlation length, ^sg and 
^cg, are obtained. This is sometimes called the dynamic 
correlation length. 

We prepare a sufficiently large lattice and calculate the 
correlation functions with high precisions, typically up to 
an order of 10^^. The correlation functions asymptoti- 
cally exhibit a single exponential decay when r is large 
enough. We estimate the correlation length using data 
in this asymptotic region. Figure^shows the lin-log plot 
of fsgir) and fcgif) in the Ising and the Heisenberg mod- 
els near the critical temperature. Linearity is good when 
the distance is larger than a typical distance, r„iin. This 
behavior guarantees our single-exponential fitting when 
estimating the correlation lengths. Values of Tmin are 
summarized in Sec. IIII El 

Correlation functions sometimes exhibit bending to 
convergence as shown in Fig.QJb). It is considered to be 
an outcome of the boundary effect. Since we work with 
the skewed periodic boundary conditions, the farthest 
distance is L/2. The correlation data near r = L/2 = 29 
are influenced by the boundary conditions. Therefore, 
the spin-glass data when r > 25 show convergence. We 
cut off such data when we estimate the correlation length. 
We also cut off data which fluctuate from the fitting line. 
The minimum distance and the cut-off distance are de- 
termined so that the relative deviations of data from the 
fitting line stay within 5%. Arrows in Fig. ^ depict the 
cut-off distance. 

It is also noted that the chiral-glass correlation func- 
tions are smaller than the spin-glass ones by three digits. 
Estimations of the chiral-glass correlation length are usu- 
ally very difficult, and the error bars become relatively 
large. 



III. METHOD 

A. Nonequilibrium relaxation method 

We start the simulations with random spin configura- 
tions, which is a paramagnetic state at T = oo. The 
temperature is quenched to a finite value at the first 
Monte Carlo step. We calculate the physical observ- 
ables at each observation time, and obtain the relaxation 
functions. Changing the initial spin state, the random 
bond configuration, and the random number sequence, 
we start another simulation to obtain another set of the 
relaxation functions. We calculate the average of the re- 
laxation functions over these independent Monte Carlo 
simulations. 

It is necessary to check a time scale when the finite- 
size effects appear in the relaxation functions. We must 
stop the simulations before the size effects appear. The 
nonequilibrium relaxation method is based upon tak- 
ing the infinite-size limit before the equilibrium limit is 
taken. Data presented in this paper can be regarded as 
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FIG. 1: (Color online) Spin- and chiral-glass correlation func- 
tions are plotted versus distance r at various time step t. 
(a) Results of the Ising model. The temperature T — 1.2 is 
near the critical point. The lattice size L = 71 for t = 10"^ 
and L = 87 for t > 10'*. (b) Results of the Heisenberg 
model. The temperature T = 0.21 is near the chiral-glass 
transition temperature. The lattice size L = 59. Arrows de- 
pict a cut-off distance when estimating the correlation length. 

those of the infinite-size system. 



B. Finite-time scaling analysis 

We estimate the transition temperature using the 
finite-time-scaling analysisi^ This is a direct interpreta- 
tion of the conventional finite-size-scaling analysis consid- 
ering the dynamic scaling hypothesis: t ~ Through 
this relation, the "size" is replaced by the "time." Here, 
the temperature dependence of the noncquilibrium dy- 
namic exponent, Zs{T), written in Eq. ^ should be taken 
into account. The temperature coefficient, l/zsgTgg, is 
obtained beforehand by the scaling analysis of the dy- 
namic correlation length. The following equation repre- 
sents the scaling formula for the spin-glass susceptibility. 

Xsg(t)i-^-/^=(^'''- = xss{t/{T - r,g)-^=(^)''-) (15) 



Here, x denotes a universal scaling function. We 
plot data of the spin-glass susceptibility multiplied by 
f-j,Jz,(T)u,^ against i/(r-rsg)~''='^^^ We choose scal- 
ing parameters, Tgg, 7sg, and Vsg so that they yield the 
best universal scaling function. 

We perform a scaling analysis using a criterion which 
we proposed in the previous paper^ In general, it is diffi- 
cult to choose proper scaling parameters because the uni- 
versal scaling function form is not known. One usually 
approximate the scaled data by some polynomial func- 
tions, and choose scaling parameters so that the fitting 
error becomes smallest. However, if we change the scal- 
ing parameter, the function form itself changes. Then, a 
significance of the fitting error changes. A value of the 
fitting error sometimes misleads us to incorrect results. 
Therefore, we take another strategy. 

Since the correct scaling function is universal, it is free 
from human factors. It should be robust against a change 
of a temperature range of data we use in the analysis; it 
should be robust against a change of discarding steps in 
simulations; it should be robust against a change of pro- 
cedures determining scaling parameters. In this paper, 
we consider the following three procedures: 

• (Proccdurc-1) For a given value of 7sg/^'sg, search 
Tgg and I'sg that yield the smallest fitting error. 

• (Procedure-2) For a given value of 7sg, search Tgg 
and Vsg that yield the smallest fitting error. 

• (Proccdurc-3) For a given value of i^sg, search Tgg 
and 7sg that yield the smallest fitting error. 

Here, we mean the smallest fitting error by the least- 
square deviations when we approximate the scaled data 
with a polynomial function of seven orders. We also 
change the temperature range and a number of discarding 
steps for each procedure. Scaling results systematically 
depend on these procedures. The most probable result 
is obtained by a crossing intersection or a merging point 
because it should be independent of the procedures. An 
error bar of the result is given by a range of scaling pa- 
rameters which yield a good scaling plot. 

C. Replica number VS Sample number 

A replica number controls an accuracy of the thermal 
average. A sample number controls an accuracy of the 
configurational average. In order to improve an accuracy 
of the thermal average, m should be large. However, if we 
set it large, a sample number is restricted, which causes a 
poor configurational average. So, there is a question. In 
order to improve an accuracy of the final results, which 
number should be set large first, a replica number m or 
a sample number iVg? Most of the simulational studies 
in spin glasses take a choice of m = 2. 

Figure |21 shows the answer. Relaxation functions of 
the spin-glass susceptibility in the Ising model are plot- 
ted both in the paramagnetic phase (a) and in the spin- 
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FIG. 2: (Color online) Relaxation functions of the spin-glass 
susceptibility of the Ising model when a replica number m and 
a bond sample number A's are changed keeping the product 
almost constant, m x A^s ~ 300. (a) T = 1.35 > Tc, and 
(b) T = 0.75 < Tc- Reference curves are plotted with black 
dotted line, which are obtained by (a) L = 83, m = 31, and 
iVs = 86; (b) L = 71, m = 128, and = 198. Each function 
is shifted to the t direction in order to compare the deviations. 



glass phase(b). A product m x A^g is set almost equal to 
300 and the lattice size is set to L = 19. A total num- 
ber of spin updates and a total computational time are 
same. We also plot a reference curve which is obtained 
by a sufficiently large simulation. All relaxation curves 
of L = 19 are consistent with the reference curves. Sta- 
tistical fluctuations are large when the replica number 
is small. They are suppressed as we increase a replica 
number and decrease a sample number. However, when 
a sample number is one, there is a systematic deviation 
from the reference curve. This deviation depends on each 
random-bond sample. The dependence becomes small 
and the deviation appears later if we increase the lattice 
size. A sample number can be set small in the nonequi- 
librium process of a sufficiently large system, where the 
size can be considered as infinity. Therefore, a replica 
number should be set larger first when the total com- 
putational time is restricted. We need at least several 
sample numbers to avoid systematic deviation due to the 
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FIG. 3: (Color online) A replica overlap in the Heisenberg 
model between the exactly same replica with the same ran- 
dom number sequence but with a different usage of a square- 
root function. Replica (a) uses a built-in sqrt function, and 
replica (b) uses an exponential function. A heat-bath algo- 
rithm and a Metropolis algorithm are compared. A temper- 
ature, T/J — 0.21, is near the chiral-glass transition point. 



sample dependence. The statement is true both in the 
paramagnetic phase and in the spin-glass phase. 

In this study, we set m = 128 and m = 192 for the 
Ising model, to = 32 and rn ~ 64 for the XY model, and 
TO = 32 for the Heisenberg model. Simulation conditions 
are summarized in Sec. IIII El 



D. Updating probability: Metropolis VS Heat-bath 

Spin update dynamics in this paper is the Metropolis 
algorithm. A heat-bath algorithmic is usually adopted 
in a simulation of the Heisenberg model. However, we 
found a chaotic behavior in a noncquilibrium dynamics 
of the heat-bath algorithm as shown in Fig. |21 

We prepare two real replicas, (a) and (b), with the 
same bond configurations with the same initial spin 
configurations with the same random number sequence. 
Only one difference is a use of a square-root function. A 
replica-(a) uses a built-in sqrt(x) function for ^/x, while 
a replica-(b) uses an exponential function x**0.5. An 
spin overlap function between two replicas is observed. 
Mathematically, it must stay unity forever. However, the 
overlap decays after some characteristic time in the heat- 
bath algorithm. This time depends on the temperature 
and also a compiler option which controls numerical pre- 
cisions. It becomes longer for higher precisions. To the 
contrary, an overlap in the Metropolis algorithm stays 
unity. 

Mathematical functions are numerically approximated 
in computer calculations. A square-root function, ^/x, is 
one of the most difficult function to approximate. The 
error is negligible for each function call. However, suc- 
cessive calls may enlarge the error. In the heat-bath 
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algorithm!^ a spin state, Si, is a function of a local 
molecular field from the neighboring spins, Sj. The 
neighboring spin state, Sj , is also a function of the molec- 
ular field which contains Si. Therefore, a spin state Si 
given by a complex function of spin states of itself and 
neighboring spins in the past. A square- root function is 
called successively in this procedure. Therefore, a small 
numerical error may be enhanced. A spin state may be- 
come different from the correct one. Wc cannot know 
which state is correct because every computer simulation 
contains a numerical error. 

In the Metropolis algorithm, each spin state is gen- 
erated independently by random numbers. There is no 
correlation between spin-state candidates in a time se- 
quence. The overlap in Fig. stays unity as mathemat- 
ically required. In this paper, wc study a spin-chirality 
separation in the Hciscnberg model. The separation may 
be observed in a very narrow temperature range. We deal 
with this very naive problem by the relaxation method. 
Therefore, we use a Metropolis algorithm in order to ex- 
clude out the unphysical chaotic behavior from the relax- 
ation process. 



E. Simulation conditions 

In this subsection, we summarize simulation conditions 
in Tabs, n ni Uni HVl and |Vl A sample number N^, a 
linear size L, a minimum distance of the spin correla- 
tions Tmin vary with a Monte Carlo time step. Numbers 
above right arrows in tables denote a time step when the 
condition changes. 

In the Ising model, a replica number is 128 except for a 
case of T = 1.4 and L = 111, where it is 192. Simulation 
conditions when the dynamic correlation length is esti- 
mated are listed in Tab. and those when the spin-glass 
susceptibility is calculated are listed in Tab. ^1 

In the model, a replica number is 64 when the dy- 
namic correlation length is estimated, and it is 32 when 
the susceptibility is estimated. The minimum distance 

8 for 



when estimating the correlation length is 
the spin glass and rjjfjjj = 5 for the chiral glass. Sim- 
ulation conditions when the dynamic correlation length 
is estimated are listed in Tab. IIIII and those when the 
spin-glass susceptibility is calculated are listed in Tab. 

In the Hciscnberg model, a replica number is 32. Sim- 
ulation conditions arc hsted in Tab. El 



IV. RESULTS OF THE ISING MODEL 

A. Dynamic correlation length 

Relaxation functions of the spin-glass correlation 
length are shown in Fig. 0] Algebraic behavior is ob- 
served in the nonequilibrium relaxation process at all 
temperatures. A nonequilibrium dynamic exponent is 
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TABLE I: Simulation conditions of the Ising model 
when is estimated. Numbers above right ar- 

rows denote a time step when the value changes. 
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TABLE II: Simulation conditions of the Ising model 
when Xsg is calculated. Numbers above right ar- 
rows denote a time step when the value changes. 
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TABLE III: Simulation conditions of the XY model 
when Csg/cg is estimated. Numbers above right ar- 
rows denote a time step when the value changes. 



evaluated from a slope of the relaxation function. It is 
well approximated by 



with ^0 = 0.54, r = 4, and 

1.18 



ZsiT) = (5.1 ±0.1) X 



T 



(16) 



(17) 



A product, Zsg^lsg, is a fitting parameter in this expres- 
sion. It is estimated by the following scaling analysis. 
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TABLE IV: Simulation conditions of the XY model 
when Xsg/cg are calculated. Numbers above right ar- 
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TABLE V: Simulation conditions of the Heisen- 
berg model. Numbers above right arrows de- 

note a time step when the value changes. 



The relaxation functions of the correlation length 
are scaled to a single line if we plot them against 
T/Tgg ln(t/T). This plot is proposed by Bert et ali for 
the experimental analysis. A product Zsg^sg is obtained 
first from a slope of the universal scaling line when we 
plot the data against Tln(</T). Here, only a character- 
istic time r is the scaling parameter, which is very small 
in the Ising model: t = 4. Using a value of ZggTsg we 
perform the finitc-time-scaling analysis in the next sub- 
section. Then a value of Tgg = 1.18 is obtained. Us- 
ing the estimated critical temperature we obtain the dy- 
namic critical exponent at the transition temperature, 
Zsg^ = (5.1 ± 0.1). The result is shown in Fig. El The re- 
laxation functions in the paramagnetic phase show a con- 
vergence crossover. It appears later as the temperature 
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FIG. 4: (Color online) Relaxation functions of the 
spin-glass correlation length in the Ising model. The 
temperatures are denoted aside the fitting lines. 
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FIG. 5: (Color online) The spin-glass correla- 

tion length in the Ising model is plotted against 
r/rsgln(f/r) using r = 4 and T^g = 1.18. A fit- 
ting line with Zsg — 5.1 is plotted by dotted line. 



approaches Tgg. Wc only consider the nonequilibrium 
process before the crossover time when evaluating the 
nonequilibrium dynamic exponent. A temperature de- 
pendence of Zs (T) is smooth above and below the critical 
temperature. This plot suggests that the nonequilibrium 
process is independent of the critical temperature. 



B. Finite-time-scaling analysis of Xsg 

The relaxation functions of the spin-glass susceptibility 
Xsgit) are plotted in Fig. Wc perform the finite-time 
scaling analysis using these data. We discard data in 
the initial relaxation regime before the tdisc step, which 
are depicted in Fig. A typical value for i^isc varies 
between 100 and 1000. The initial relaxation is very short 
compared to the equilibrium relaxation. 
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FIG. 6: (Color online) Relaxation functions of the spin- 
glass susceptibility in the Ising model. Relaxation data 
after fdisc (red curves) are used in the scaling analysis. 
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FIG. 7: (Color online) Estimates of the finite-time- 
scaling parameters in the Ising model. (a) Tsg ver- 
sus 2 — 77sg. (b) Tsg versus v^^. For each tempera- 
ture range, we plot results of two different procedures 
(circle and square symbols) and two different discarding 
steps (filled and unfilled symbols). Rectangles depict error 
ranges, which covers most of possible scaling parameter sets. 
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FIG. 8: (Color online) Finite-time-scaling plots in 
the Ising model. (a) A scaling plot of Xsg us- 
ing the obtained results. (b) A scaling plot of 
the correlation length in the equilibrium convergence. 



We follow the procedures explained in Sec. lIII Bl Three 
different temperature ranges are considered in this anal- 
ysis. One consists of all the 17 temperatures of Fig. 
ranging 1.3 < T < 1.6. The second one consists of 
the lower 12 temperatures ranging 1.3 < T < 1.46. 
The third one consists of the lower 8 temperatures rang- 
ing 1.3 < T < 1.42. For each temperature range, we 
adopt three procedures of estimating the scaling param- 
eters. We find that the procedure-2 and the procedure-3 
yield consistent results. Therefore, we compare results 
of procedure-1 and procedure-2. Wc also change the dis- 
carding time step, which is distinguished in Fig. [T] with 
filled and unfilled symbols. Unfilled symbols depict data 
when the first tdisc steps are discarded, where tjisc are as 
depicted in Fig. Filled symbols depict data of discard- 
ing the first 3tdisc steps. 

Figure [Tfa) shows the obtained critical temperature 
plotted versus 2 — r]sg. Behaviors of the estimated crit- 
ical temperature depends on the procedures, the tem- 
perature ranges, and the discarding steps. These data 
seem to cross or merge at the most probable parameter 
point of (2 — 77sg,Tsg) ~ (1.96, 1.18). We put error bars 
so that they cover most of possible scaling parameters. 
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Figure Cfb) shows the critical temperature plotted ver- 
sus the critical exponent z/gg. The scaling results cross or 
merge at {vsg,Tsg) ~ (1-73, 1.18). Two estimates for the 
transition temperature are consistent. An exponent 7sg 
is estimated from Tjsg and Vsg- Using the estimated scal- 
ing parameters we plot the scaling functions in Fig.|Hta). 
We also plot in Fig.lH^b) a result of the finite-time scaling 
analysis of the dynamic correlation length. It yields the 
following scaling relation: 



-l/^s(T) _ 



e;g(t/(T-T.g)-^=(^)'^-) 



(18) 



The shape of the universal scaling function is robust 
against changes of discarding steps and the temperature 
ranges. 

From the present scaling analysis, we conclude the crit- 
ical temperature and the critical exponents of the Ising 
model in three dimensions: 



Isg 
^sg 
Vsg 
Zsg 



1.18 ±0.02 
3.4 ±0.2 
1.75 ±0.10 
0.04 ±0.03 
5.1 ±0.1. 



(19) 
(20) 
(21) 
(22) 
(23) 



Previous estimates^-^^i^S-^Si^ii^ for T^g he between 1.1 
and 1.2. They are roughly categorized into two groups. 
Group I: Tgg close to 1.1 and i^sg close to 2m^i^ Group II: 
Tgg close to 1.2 and i^sg close to 1.342il^i^ Our estimate 
for Tsg belongs to the latter group. On the other hand, 
an estimate for Vsg is rather close to the former group and 
agrees with the experimental result 4^ which is i^gg ~ 1.7. 
It is very hard to determine a value of Vsg. 

Recently, Campbell et. alm^ claimed that a choice of 
a scaling variable may affect the scaling results. Partic- 
ularly, a value of ly changes. We applied their method 
and checked the claim. A scaling plot is possible with 
Tsg = 1.18, i^sg — 2.18, and 7sg ~ 4.15, when we use a 
scaling vaiable {Psg — P) instead of (T — Tsg). The crit- 
ical temperature and the anomalous exponenet Tysg are 
not affected by this change. We consider that these two 
estimates are reliable. To the contrary, two exponents, 
Vsg and 7sg, become larger than our original estimates 
and are consistent with the values of group I. We cannot 
exclude out either estimate within our amount of simu- 
lations. 



V. RESULTS OF THE XY MODEL 

A. Brief introduction to the vector models 

It is widely accepted that the spin-glass transition oc- 
curs in the Ising model i^24^iiSi^SiiLi£ However, spins 
of many spin-glass materials are of the Heisenberg 
type. Early numerical investigationa^SiSiiS concluded 
that there is no spin-glass transition in this model. 
Kawamurai^SiSLSSi^ proposed the chirality-mcchanism in 



order to explain the real spin-glass transition. The spin- 
glass transitions in real materials occur driven by the 
chiral-glass transition. This scenario assumes that there 
is no spin-glass transition in the isotropic model42*^ Re- 
cently, there are several investigations which suggest the 
simultaneous transitioni^iiSi^liSi^ 

The situation is same in the XY spin-glass model 
in three dimensions. A domain-wall energy analysis^ 
and a Monte Garlo simulation analysis^ found no 
spin-glass transition. However, a possibility of a fi- 
nite spin-glass transition has been suggested by recent 
investigation82li^SiS24^2iS£. Since the chirality is defined 
by spin variables, the chiral-glass transition occurs triv- 
ially if the spin-glass transition occurs. Therefore, it is 
crucial to check whether or not both transitions occur 
at the same temperature. This is a very naive prob- 
lem. For example, the spin-chirality separation is ob- 
served in the fully frustrated nonrandom XY model in 
two dimensions However, the difference between the 
critical temperatures is only 1%. 



B. Dynamic correlation length 

Relaxation functions of the spin- and the chiral-glass 
correlation length are shown in Figs. a) and (b). Al- 
gebraic behaviors arc observed in the nonequilibrium re- 
laxation process both above and below the critical tem- 
perature. Slope of relaxation functions of the spin-glass 
correlation length continuously depends on the tempera- 
ture. There is no anomaly at the transition temperature, 
Tsg = 0.435, which will be estimated in the next subsec- 
tion. This is the same behavior with the Ising model. On 
the other hand, relaxation functions of the chiral-glass 
correlation length do not depend on the temperature 
at higher temperatures (T > 0.45). It changes the be- 
havior at T ~ 0.45 becoming temperature-dependent at 
lower temperatures. This is a distinct difference between 
the spin- and the chiral-glass nonequilibrium dynamics. 
From the scaling analysis of the correlation lengths, we 
obtain the temperature dependences of the nonequilib- 
rium dynamic exponents. Figurc|HIc) shows results of the 
scaling. The scaling is possible both below and above Tsg 
for the spin-glass correlation length, but it is only possi- 
ble below Tsg ~ Teg = 0.45 for the chiral-glass one. It is 
also noted that a characteristic time is common between 
the spin and the chirality, and it is larger than the Ising 
one. 

We obtain temperature dependences of the dynamic 
exponent for the spin glass (2:s(T)) and the chiral glass 
(zc(T)) as follow. 



^c(T) = 
^ciT) = 

Here, we use Ts< 



(5.0 ±0.1) X 
(5.1 ±0.2) 
(5.1 ±0.2) X 
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T 
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(T > 0.45) (25) 
(T < 0.45) (26) 
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FIG. 9: (Color online) Relaxation functions of the 
correlation lengths in the XY model. (a) The spin- 
glass correlation length. (b) The chiral-glass correla- 
tion length. (c) Scaling plots of ^sg and ^cg versus 
r/rcln(t/r). Scaling parameters are as denoted in the fig- 
ure. The chiral-glass scaling is only possible at T < 0.45. 
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FIG. 10: (Color online) Relaxation functions of the 
spin- and the chiral-glass susceptibility in the XY model. 



when the spin-glass exponent equals to the the chiral- 
glass one. The spin-glass transition occurs at this tem- 
perature. Both exponents behave in the same manner 
below the temperature. 



Finite-time-scaling analysis of Xsg and Xc 



Figure 1101 shows the relaxation functions of the spin- 
and the chiral-glass susceptibility. Some of the data are 
taken from our previous worki^ We added several tem- 
perature points and increased a total Monte Carlo step 
particularly near the transition temperature. In the pre- 
vious work, we supposed that the dynamic exponent is 
independent from the temperature. Crossing behaviors 
of the scaling results were not good. Consequently, the 
error bars of obtained results became large. These fail- 
ures are possibly due to the temperature dependence of 
the nonequilibrium dynamic exponent. In the present 
analysis we use Eqs. (|24|l and l|25|) . A scaling procedure 
is same as the Ising case. 

Figures [TTl (a) and (b) show estimations of the scaling 
parameters. We determine the best estimates from the 
figure so that the results are robust against the difference 
of procedures. For the spin-glass transition: 



estimated in the next subsection. A coefficient to the 
temperature-ratio term is the dynamic critical exponent 
at the transition temperature. It agrees with each other 
within error bars. Therefore, the dynamic critical expo- 
nents of the spin-glass and the chiral-glass transitions are 
considered to take the same value. The nonequilibrium 
dynamic exponent of the chiral glass changes behavior 



T — 

J^sg = 

»7sg = 

2sg = 



0.435 ±0.015 
2.30 ±0.25 
1.25 ±0.15 
0.18 ±0.07 
5.0±0.1. 



(27) 
(28) 
(29) 
(30) 
(31) 
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FIG. 11: (Color online) Estimations of the finite-time- 
scaling parameters in the XY model, (a) Tc versus 2 — 7]. 
(b) Tc versus v. Spin-glass (chiral-glass) results are plot- 
ted against the lower(upper) x-axis. The circle symbols 
are results of the scaling procedure- 1, and the square sym- 
bols are those of the procedure-2. Unfilled symbols depict 
data when the first 400 steps are discarded. Filled sym- 
bols depict data when the first 1000 steps are discarded. 

For the chiral-glass transition: 

Teg = 0.445 ±0.015 (32) 

7cg = 1.65 ±0.20 (33) 

ly^s = 1.05 ±0.10 (34) 

77cg = 0.40 ±0.05 (35) 

2cg = 5.1 ±0.2. (36) 

The spin- and the chiral-glass transition temperatures 
cannot be distinguished to be different. The chiral-glass 
transition occurs at a little higher temperature but the 
difference is within the error bars. A value of the dynamic 
critical exponent coincides with each other within the er- 
ror bar. Both transitions are considered to be dynami- 
cally universal. On the other hand, the static exponents 
take different values. 

The scaling plot using the obtained results is shown in 
Figs. I12r a) and (b). All data suitably fall onto a sin- 
gle scaling function. The obtained transition tempera- 



FIG. 12: (Color online) Finite-time-scaling 

plots of (a) Xsg and (b) Xcg using the esti- 
mated scaling parameters in the XY model. 



ture is between our previous estimate^* supposing a con- 
stant dynamic exponent and other estimatesi^SiSS Our 
previous estimates are Tgg = 0.455(15), 7sg = 1-7(2), 
^^sgt'sg = 4.8(4) for the spin-glass transition and Teg = 
0.467(10), 7cg = 1.4(1), and Zcgi'cg = 4.7(2) for the 
chiral-glass transition. The results of the spin-glass tran- 
sition are strongly influenced by the temperature depen- 
dence of Zs (T) . Granato^S provided the following for the 
spin-glass transition: Tgg = 0.39(2), v = 1.2(2), and 
z = 4.4(3). Kawamura and Li^ provided the following 
for the chiral-glass transition: Teg = 0.39(3), v = 1.2(2), 
•y/iy = 1.85(20), and z 7.4(10). These transition tem- 
peratures are lower then our present estimates. 



VI. RESULTS OF THE HEISENBERG MODEL 

A. Dynamic correlation length 

Relaxation functions of the spin- and the chiral-glass 
correlation lengths are shown in Figs. Iiat a) and (b). We 
observe algebraic increases of both correlation lengths. 
We estimate the dynamic exponent using these data. 
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FIG. 13: (Color online) Relaxation functions of the cor- 
relation lengths in the Heisenberg model, (a) The spin- 
glass correlation length. (b) The chiral-glass correla- 
tion length. Fitting lines are plotted with broken lines, 
(c) The scaling plot of plotted versus T/Tsg ln(f/r). 



The chiral-glass exponent may have a rather large er- 
ror bar because numerical fluctuations of the chiral-glass 
data are large. Evaluatios of the chiral-glass quantities 
are generally difficult because the chiral-glass correlations 
are much smaller than the spin-glass ones as shown in 
Fig. n^b). It is necessary to increase a sample number 
in order to suppress the fluctuations. Within the present 
computational facilities, it is not realistic to increase a 



Monte Carlo step number until the chiral-glass correla- 
tion length grows comparable to the present spin-glass 
data. For example, we need a hundred times longer steps 
to wait until ^cg reaches five lattice spacings. 

It is also noted that the algebraic increase in the low- 
temperature phase appears later as the temperature de- 
crease. It is considered to be an outcome of the slow 
dynamics of the spin-glass phase. A time range to es- 
timate the dynamic exponent is restricted to be short. 
Reliability of the estimated value may become poor. 

Even though the numerical fluctuations are large, we 
observe a temperature-dependence of the dynamic ex- 
ponent as in the Ising model and the XY model. Loga- 
rithmic slope of the spin-glass correlation length depends 
on the temperature irrespective of the transition tem- 
perature. That of the chiral-glass correlation length de- 
pends on the temperature only when T < 0.18. These 
evidences are checked by a scaling plot of Fig. I13r c). 
The spin-glass data can be scaled to a single line both 
above and below the critical temperature, Tgg = 0.18, 
which will be estimated in the next subsection. On 
the other hand, the chiral-glass data are scaled only be- 
low Tsg. The high temperature data systematically de- 
viate from the scaling line. They are better approxi- 
mated by a temperature-independent algebraic function: 
^cg = 0.62 X (i/3000)i/^ ". Therefore, the chiral-glass 
dynamic exponent is constant above Tgg and becomes 
temperature dependent below Tsg- The temperature de- 
pendence of the chiral-glass dynamic exponent below Tgg 
is same as the spin-glass one. These behaviors are same 
as the XY model. A difference is a characteristic time r, 
which is much larger than the Ising and the XY models. 

The temperature dependences of the nonequilibrium 
dynamic exponent in the Heisenberg model, Zs{T) for 
the spin glass Zc(T) for the chiral glass, are summarized 
as follow. 



(5.1 ±0.1) X 
(5.0 ±0.2) 
(5.1 ±0.2) X 



0.18 



0.18 



(37) 

{T > 0.18) (38) 
(T < 0.18) (39) 



Values of the nonequilibrium dynamic exponent at T = 
Tgg = 0.18 agree with each other. This is the temper- 
ature where the chiral-glass dynamic exponent changes 
the behavior. It can be considered that dominant dy- 
namics changes when the T-dependent Zs{T) catches up 
with the T-independent Zc{T) at T = 0.18. The spin- 
glass nonequilibrium dynamic exponent, Zs(T), is larger 
than the chiral-glass dynamic critical exponent, Zcg, be- 
low the temperature. Then, the dynamics of the chirality 
obeys that of the spin because the chirality is defined by 
the spin. As the result, Zc{T) behaves in a same manner 
with Zs{T). 
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FIG. 14: (Color online) Relaxation functions of 

the spin-glass and the chiral-glass susceptibility in 
the Heisenberg model. Data of the chiral-glass sus- 
ceptibility are multiplied by 535 in order to show 
in the same scale with the spin-glass susceptibility. 



B. Finite-time-scaling analysis of Xsg and Xcg 



Figure im shows relaxation functions of the spin- and 
the chiral-glass susceptibility. Data of the chiral-glass 
susceptibility are multiplied by 535 in order to show them 
in a same scale with the spin-glass susceptibility. Am- 
plitude of the chiral-glass susceptibility is roughly 1000 
times smaller than the spin-glass one. The chiral-glass 
order is very small. 

Figures El(a) and (b) show estimations of scaling pa- 
rameters. We obtain the transition temperature and the 
critical exponents, where the scaling results are robust 
against the change of procedures. It is clear from the 
figures that the spin-glass transition and the chiral-glass 
transition occur at different temperatures. Error bars of 
the critical exponents are rather large compared to those 
of transition temperatures. Particularly, a value of Vsg 
systematically shifts to a higher side when we decrease 
the lowest temperature of data used in the scaling anal- 
ysis from T = 0.215 to T = 0.1925. Then, a value of ry^g 
is influenced. However, a value of 7sg does not change. 
Our estimates for the spin-glass transition are 



7sg 

Vsg 

Zap- 



0.18 ±0.01 
2.9 ±0.2 
1.55 ±0.15 
0.1 ±0.1 
5.1 ±0.1, 



(40) 
(41) 
(42) 
(43) 
(44) 



0.24 
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^0.21 
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0.16 
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FIG. 15: (Color online) Estimations of finite-time-scaling pa- 
rameters in the Heisenberg model, (a) Tsg/cg versus 2 — 77ag/t.g. 
(b) Tsg/cg versus «^sg/cg- Spin-glass(chiral-glass) results are 
plotted against the lower(upper) axis. The circle symbols 
are results of the scaling procedure- 1, the squares are of the 
procedure-2, and the triangles are of the procedure-3. Un- 
filled symbols depict data of 1000 discarding steps, and filled 
symbols depict those of 3000 steps. Color distinguishes the 
temperature range of data used in the scaling analysis. Re- 
sults of 0.1925 < T < 0.27 are obtained by discarding 3000 
steps (unfilled symbols) and 10000 steps (filled symbols). 



and those for the chiral-glass transition are 



7cg 

J^cg 



0.23 ±0.01 
0.65 ±0.15 
0.70 ±0.15 
1.1 ±0.1 
5.0 ±0.2. 



(45) 
(46) 
(47) 
(48) 
(49) 



In our previous work^^i we assumed that the dynamic 
exponent does not depend on the temperature. The ob- 



tained results are Ts^ = 0.21 



+0.01 
-0.03 



and Tc, 



0.22 



+0.01 
-0.04- 



Present estimates roughly correspond to the upper and 
the lower bound of the previous estimates. The error 
bars shrink, and two transitions are resolved to be differ- 
ent, which is made possible taking into acount the tem- 
perature dependence of the spin-glass dynamic exponent. 
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FIG. 16: (Color online) (a)A finite-time-scaling plots 
of Xsg using the estimated parameters. Two possible 
chioces of scaling parameters are shown. As we dis- 
card more steps, a value of v^^ increases. (b) A scal- 
ing plot of Xcg- Data of the chiral-glass phase, Tsg < 
Teg, are also scaled using the same parameters. 



Therefore, the spin-glass results are much influenced. Re- 
sults of the chiral-glass transition are consistent with the 
previous ones, because an assumption of a constant dy- 
namic exponent is good in this case. Even though the 
static exponent, Vsg and 7sg, take different values, Zgg and 
Zcg coincide within the error bars. They also agree with 
the dynamic critical exponents of the Ising and the XY 
models. The spin-glass and the chiral-glass transitions 
in three dimensions can be considered to be dynamically 
universal. 

Scaling plots using the obtained results arc shown in 
Figs. llGf a) and (b). For the spin-glass transition, we 
plot two results with different choices of scaling parame- 
ters. They are within the error bars of our estimates. A 
value of Vsg tends to increase when the lowest tempera- 
ture of the scaling data decreases and when we increase 
the discarding steps. The transition temperature and a 
value of 7sg are rather robust against these changes. A 
scaling plot using a parameter set with = 1.48 shows 
some systematic deviations in the long time limit at low 
temperatures. Therefore, our best choice in the present 



analysis is 7sg = 2.9, = 1.6, and Tgg = 0.185. A ratio 
of scaling parameter, 7sg/i'sg, agrees with the logarithmic 
slope of the spin-glass susceptibility of Fig.^^at T = Tgg. 

For the chiral-glass transition, we plot data both above 
and below the chiral-glass transition temperature. A 
scaling to a single curve is good above the chiral-glass 
transition temperature. A scaling is also possible for the 
data in the chiral-glass phase: Tgg < T < Teg. The scal- 
ing function is bending upward, which suggests that a 
ratio of critical exponent, "f/v, is increasing as the tem- 
perature decreases. 

We also applied the /3-scaling proposed by Campbell 
et. alm^ and checked our estimates. As in the Ising 
model, the critical temperature and the anomalous ex- 
ponent 77sg do not change, while Vsg and 7sg increase. 
The estimates are Tgg = 0.18, J^sg = 3.7, and 7sg = 6.4. 
Two exponents, Vsg and 7sg, become particlarly large. 
The present numerical data may be insufflcient to deter- 
mine these critical exponents. On the other hand, the 
spin-glass transition temperature is consistent with our 
original estimate. We obtained this value by two scaling 
analyses and by a temperature dependence of the chiral- 
glass dynamic exponent. Wc consider that the value is 
very reliable. 



VII. SUMMARY AND DISCUSSION 

The nonequilibrium dynamics and the spin-glass tran- 
sitions in three dimensions are clarified in detail. It is 
confirmed that the nonequilibrium dynamic exponent of 
the spin glasses continuously depends on the temper- 
ature without anomaly at the transition temperature. 
The spin-glass dynamics in the high-temperature phase 
smoothly continues to the low-temperature phase. The 
frozen-spin domain grows in the algebraic law controlled 
by Zs{T) irrespective of the critical temperature. The 
growth in the high-temperature phase stops when the 
domain size reaches the equilibrium value. The smooth 
temperature dependence of the nonequilibrium dynamic 
exponent is considered to be a common characteristic be- 
havior of the spin glasses. We will need to explain the 
origin regardless the spin-glass transition, which is left 
for the future study. 

In the vector spin models, the nonequilibrium dynamic 
exponent of the chiral glass shows a different behavior. 
It does not depend on the temperature above the "spin- 
glass transition temperature" . It remains constant until 
it changes the behavior at T = Tgg. When T < Tgg, 
the chiral-glass dynamic exponent is same with the spin- 
glass one. Therefore, the dynamics of the chiral glass 
is essentially same with the spin-glass dynamics in the 
low-temperature phase. They are different only in the 
high-temperature phase. 

Our observation suggests the following scenario. Dy- 
namics of the spin and the chirality are separated at 
temperatures above Tsg- The dynamic exponent of the 
spin glass is smaller than the chiral-glass one. The 
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spin-glass state evolves with time faster than the chiral- 
glass state. Therefore, the chiral glass freezes first at a 
higher temperature. This is the chiral-glass transition 
as Kawamurai2Si^L22i^ proposed. Since the spin and the 
chirality are separated, the spin-glass dynamic exponent 
shows no anomaly at T = Teg. The situation changes at 
T = Tgg, where the spin-glass dynamic exponent over- 
comes the chiral-glass dynamic critical exponent, Zcg. 
The spin-glass transition occurs at this temperature. The 
value, Zgg = ^cg — 5.1, is found to be common to all the 
spin-glass models in three dimensions. Below the spin- 
glass transition temperature, the spin-glass dynamics is 
slower than the chiral-glass dynamics above Tgg, because 
Zs{T) > Zcg. Since the chirahty is defined by the spin 
variables, the chirality obeys the spin dynamics once the 
spin freezes and the spin-glass dynamics becomes dom- 
inant. The spin and the chirality is coupled below Tgg, 
and both dynamic exponents show the same temperature 
dependence. Therefore, this is the spin-glass phase. 

Using the obtained temperature dependences of the 
dynamic exponent wc performed the finite-time-scaling 
analysis on the spin- and the chiral-glass susceptibil- 
ity. The transition temperature and the critical ex- 
ponents are estimated with high precisions. In the 
Ising model, the results arc compared with the previous 
worko^3,44^45,46,47,48 ^he experimental result .i^ We 

also checked that a change of a scaling variable may affect 
a value of Vsg- In the model, the spin- and the chiral- 
glass transitions occur at the almost same temperature. 
In the Hciscnberg model, the spin- and the chiral-glass 
transitions are resolve to be different. It is the lower and 
the upper bound of our previous estimated Therefore, it 
becomes clear that the chiral-glass transition without the 
spin-glass order occurs in the Heisenberg model, as the 
chirality-mechanism predictedi^Si2L^2i22ii2i^ The chiral- 
glass phase continues for Teg >T> Tgg. 

The chirality mechanism insists that the spin-chirality 
separation is only observed in the long time and in the 
large size limit 4i A crossover from the spin-chirality cou- 



pling regime at short scales to the decoupling regime at 
long scales is considered to be important. Discrepancy 
regarding whether the simultaneous spin-chirality tran- 
sition occurs or not is explained by the crossover. It is 
considered that the true transition will not be observed 
unless the simulation goes beyond the crossover length- 
time scale. However, wc are able to observe the spin- 
chirality separation near the chiral-glass transition tem- 
perature within a rather short time and a small length 
scale. It was made possible by clarifying the temperature 
dependence of the dynamic exponent in the nonequilib- 
rium relaxation process. The separation is observed as 
a different temperature dependence. If we neglect the 
effect of the temperature dependence, the spin-chirality 
separation will be smeared out. Then, wc arc lead to the 
simultaneous transition. 

All the dynamic critical exponents at the transition 
temperatures investigated in this paper take almost the 
same value. It may be the universal exponent controlling 
the spin-glass transitions. Therefore, we consider the dy- 
namical universality about the spin- and the chiral-glass 
transitions in three dimensions. 

The spin-glass dynamics around the transition temper- 
ature has exotic behaviors, which arc quite different from 
the regular systems. A single continuous phase seems to 
exist in regard with the dynamics. It should be explained 
in connection with the picture of the low-temperature 
phase. It will be a future subject. 
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